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By Robert F. Coleman 

SUMMARY 


A preliminary analysis has been made of the conditions of 
stability of free oscillations of a hinged rotor In hovering 
flight. The case analyzed is a rotor with hinges allowing 
freedom In flapping and lagging and having a completely reversible 
cyclic pitch-control system, so that a twisting moment on a blade 
moves the control stick without hindrance from spring or friction 
constraint. 

The principal results of this study are presented In the form 
of a stability chart In which the quantities specifying average 
coning angle, pitch setting, moments of Inertia, chardwlse end 
spenwlse mass distribution are combined Into generalized parameters 
In such a way that the stability condition can be plotted on a 
single chart. 

The results show that the stability is extremely sensitive to 
chordwise center of mass, and that forward movement of the center 
of mass increases the stability. It is also shown that the 
stability depends upon coning angle and pitch setting. For all 
examples In which the caning angle is determined by a balance 
between lift forces and Inertia forces, the condition for neutral 
stability corresponds to the same point of the chart. 

Several lines of attack for further theoretical work are 
suggested, which would extend the generality and strengthen the 
validity of the analysis. 
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IHTRODOCTION 


It Is lnqportant for the safety of helicopter flight that 

suitable precautions he taken to ovoid the occurrence of 
self -excited vibrations such as blade flutter. The understanding 
of flutter phenomena is therefore ■ important in the design of 
rotary vlng aircraft. Although helicopter blade flutter has a 
certain similarity to airplane vlng flutter, the analysis is 
ecmevhat different because of the occurrence of teams associated 
with the rotation of the blades. 

Rosenberg in reference 1 has given an analysis of 
helicopter blade flutter that consists essentially of adding 
centrifugal stiffening teams to the standard vlng flutter analysis. 
But rotation can also introduce inertia coupling between flapping, 
lagging, and feathering. These coupling teams, depending upon 
the average coning angle and pitch setting of the blades, have not 
previously been considered in flutter analysis. 

The purpose of the pore sent paper is to exhibit the simplest 
case that will point out these nev effects. This is done by 
treating the extreme case vhero inertia forces and air forces are 
dominant In comparison vith elastic forces. As this case is 
characterized by a lov flutter frequency, the phenomenon is 
conveniently inferred to as lov frequency helicopter flutter. 

In reference .1 the elastic forces are considered and certain 
ineortia coupling teams are Ignored, while in this paper these 
Inertia coupling teams are considered and the elastic teams are 
ignored . The next logical step in the solution of the flutter 
problem now seems to be to combine the tvo theories into a single 
scheme of analysis that includes both elastic effects and coning 
and pitch-angle effects. This further development, however, has 
not been Included in the present paper. 

It should be mentioned that the pz*esant analysis takes no 
account of the phenomenon of Btall flutter although pitch angle 
is also an impor*tant parameter In that problem. 


ANAIXSIS 


If the cyclic pitch-control system of a rotor is considered 
to be completely reversible, so that a twisting moment on a blade 
moves the control stick without hindrance from spring or friction 
constraint then, except for ineortia effects due to rotation, the 
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"blade may "be considered to have zero frequency In torsion. Also, 
if all the blades have the same physical characteristics, the 
stability of the rotor as a -whole can be obtained from a study of 
the equations of motion for a single blade. If further, for 
purposes of vibration analysis, the hinge lines of flapping, 
lagging and feathering are assumed to pass through a common point 
on the spindle axle, then the blade may be treated as a rigid 
body vith one point fixed. Tho well known Euler's equations of 
motion for this case are 


Afljj - (B - C)ayp g = 
Bay - (C - A )a> z a> 2 - L y 
- (A - B)ayjy « 


( 1 ) 


whore 
A, B, C 


principal momenta of inertia of blade about x, y, 
and z axes, respectively 


(Oyr, cjuy, cos components of angular velocity rosolvod along 

instantaneous directions of principal axes of blade 

Lg, Ly, L z components of external moments about the fixed point ' 


X, 


subscripts referring to principal axes of Inertia of 
blade. (See fig. 1.) * 


Die average values of aerodynamic and Inertia moments will be 
considered to be balanced by suitable constant torques applied at 
the hub by the drive shaft and, if necessary, by the control 
system. The deviations from the average will then appear in the 
equations governing the vibration. 


The total external moment 1^, Ly, L z is tho sum of the 
aerodynamic moment 1^', L y ', L z * and the hub torque Qy, Qy, Q z> 

Expressions for the components of aerodynamic moment 
Lj', ly', I> 2 *- are obtained by assuming that each airfoil section 

has an aerodynamic force normal to its Instantaneous velocity and 
equal in magnitude to the conventional expression for lift of an 
airfoil in steady motion (soe fig. 2) thus: 


n> 




a 4 wi 
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vhere 

F aerodynamic force -per unit span 

V Instantaneous velocity of section 

c chord of section 

a. Instantaneous angle of attach 

e. slope of lift curve 

o 

The effect of Induced velocity Is considered only 1 a the 
choice of numerical value for the slope of the lift curve. Drag 
f 01*068 have "been ignored In this Introductory treatment. 

As the Inertia terms In Euler* s equations are expressed In 
terms of ay., co z it Is desirable to express the aerodynamic 

moments In terms of these same variables, thus: 

v = V cos a =» rto 


V„ = V sin a = ra> 
z y 

vhere r is radial distance of blade element from hinge point, 
then 

* F cos a = gpV^ ca Q sin a cos a = gpcaQr^ayOg 
F as F sin a * ^pV®ca 0 sin^ a = -gpcaQr^cuy^ 

y 

The aerodynamic moments can then be written: 




y *scLr as 




J 0 2t«v 


ir^dr 


-F z rd r 


V -IV = 03^ r a J^ -Jpca^dr 


=Tz| -lpca 0 r3ar j> 


(3) 
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-where s Is the distance of the chordwise center of mass behind 
the aerodynamic center. 

Biller's equations now become: 



- (B - C).ayc z 

- (C - A)©^ 

- (A - B)tD x a> y • 


' a j Q>z J' -^pca o sr S dr + Q^. 

oo a> P -Ipca r^dr + Q 

72 Jo^ ° 7 

" 8 IV 


a r3dr + Q 
o z 


v. 

/ 


w 


where Q x , Qy, Q z are the components of constant torque required 

to balance average values of assumed aerodynamic and Inertia 
moments. These torque components may represent such things as 
weight moment, constant structural bending moment, and constant 
bungee force In the control system. 


In order to obtain solutions of equations (4) the deviations 

from the average values of ou , cd , co are assumed to be small 

x y z 

enough so that squares and products of deviations can be neglected. 


Equations (4) are linearized by putting 


°v =0 i + I \ 

“V- n y + “V 


f 


f 




(5) 


0) 


= 0 


+ al- 


and with the assumption that Q x , Qy, 
terms not containing eo^', cDy.', cd z ' 


J 

Q z are In equilibrium with 
the constant terms become 


< x i - V™* 


Sx 

A 


- X 2 Q A + ^y^z = 

1 fiPj - " 



( 6 ) 
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and the linear terms become 


V -Hi)«n y v iovy.o 


CD 1 - I_(0 CD * + ft CD * ) + H (ft CD + ft CD * ) » 0 ^ (7) 

y 2 z x x z 2 v z zy I 


V + +0 A ,) ' H 3 2 °yV " 0 


where 


T jC_=_J 
J 1 = A 




t B - A 
J 3 “ C 


= -J- J -1- pcapSrSdr 


H. 


*’*1* P “» 


r3<lr 


H. 


•*r< 


2 pca^^dr 


The quadratic terms are neglected. 

The linear terms govern the vibrations . Their solution Is 
of the form 


It 


V " “xo* 6 


cd * « cd * e^ 
y yo 


cd 1 = co * e^ 
z zo 


> (8) 


where X Is a root of the determlnantal equation 

di-HijOy 

-hPx * aftr 

m 

X 


X 


X + zf z 
X ?x - “fy 


(9) 
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Kile equation Is a cubic in X. When expanded in the form 
X.3 + gjfi + bX + o =» 6 


where 


a » Hgp z 


» - <*1 - HiK^A 2 - I 3 9, 2 ) + (I A - - 2Sfy) 

o = SlgCIj^ - EJ n^y(-I 3 £^ + H^) 


the well known conditions for stability are 

a > 0 
b >0 
c >0 


c < ab 

The fourth condition is the important one for oscillatory 
instability. 


For the critical condition corresponding to c = ab the 
cubic equation becomes 


(X + a)(X 2 + b) = 0 



-a 

tijb" 



( 10 ) 


Oils condition represents the border line between damped and 
self -excited vibration of the system. A chart has been devised 
to show the values of the pertinent parameters, corresponding to 
this boundary between damped and self -excited vibrations. 

Bie explicit equation corresponding to c = ab 1 b 
Oydl » %) 212 (-1^ + Kfy) 

- Eg [di - hi)(ia 2 - + (i A - W* 1 #* - ^3^2] (11) 
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The average pitch setting 9 and coning angle ft are now 
Introduced "by the substitution 


Bin p 


fly «fl cos P sin 6 7 

flg ■ fl cos P cos 6 

In. terms of p end 0 equation ( 12 ) becomes 

\ 


H, - I, 


o 

cos 6 
tan® P 


. V£ 


sin 6 
tan P 




1 - tunr 0 


(A,l 

V % i 


21, tan P 
~q ~ Hg sin 0 


& 


(13) 


(14) 


STABILITI CHART 


A chart has been plotted in which all variables of equation (14) 
have been Included subject only to the restriction 

* H 3 (15) 

which is a very good approximation for a long thin body like a 
rotor blade. In terms of variables defined by 

t - % - I t cos 2 0 

I 3 tan 2 3 

Y = sin 0 
12 ten P 
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equation ( 14 ) becomes 


I „ £ - 1 ^ gY - jj 

1 - tan 2 0 (3 --f) (16) 

The chart, figure 3 , la then a plot of Y against X for 
constant values of 0 . 


DISCUSSION OF RESULTS 


The variables X and Y of equation (16) have been defined 
In such a t ray as to make the chart nearly a universal single curve. 
A sight dependence upon Q Indicated by the different curves of 
the chart, shows that It Is, In reality, a family of curves; but 
that all the curves are close together for any reasonable range of 
values of f>. The main dependence of the flutter condition upon 
the physclal characteristics and the operating condition of the 
blr.de Is consequently Implied In the form of expression defining 
X and Y. 

The most Important blade characteristic is the chord vise 
position of center of mass with respect to aerodynamic center. 

This dietance s appears In the definition of H-j. and 

consequently In X. For typical parameters, changing the center 
of mass from 1 percent ahead to 1 percent behind the aerodynamic 
center will change the value of X from -15 to 15 if 
*1 «* 0 and from -45 to -15 if 1^ = 1. The value of I-j_ 

would be close to zero for a blade with a heavy spar of circular 
cross section. It would be close to 1 for a blade with Its 
mass well distributed In the chordwlse direction. 


The variable Y is equal to the ratio of the average aero- 
dynamic moment to the average inertia moment about, the axle 
parallel to the blade chord. 


For 


average aerodynamic moment. 



average inertia moment, 


(C 


i* _ Hp sin 9 

(C - A) I2 tan P 
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The value of Y for blades vlth a flapping hinge will therefore 
normally he close to 1, so that the critical point for nearly all 
typical applications wiUr he close to »0.- Departures 

from the value Y = 1 will he associated with such things as 
gravity moment or structural bending moment in the blade, or 
possibly transient flapping conditions as in gusts. These effects 
may thus make the blade flutter if they Increase the oonlng angle 
or make it more stable if they decrease the caning angle. It 
seems, for example, as though the gravity moment would tend to make 
the blade more stable. Cases could thus be imagined that, for a 
constant pitch setting were stable at very low rpm, where gravity 
has a large relative effect, and unstable at operating rpm, where 
a typical value of Y lo 1.04. A model tested upside down would 
be expected to show the opposite effect. 


FUMHER KEFUTEMEHTS OF THEORY 


This treatment is to be considered as only a preliminary 
study of blade flutter. Many simplifying assumptions have been 
made in order to obtain a simple stability chart and for this reason 
the reeults should be applied with caution. This analysis can now 
be used as a starting point for further refinements. Seme effects 
which might be considered are: 

1. dreg forces 

2. case where I g £ 

3. elastic and friction forces 

4. unsteady- lift functions 

5. constraints such as completely irreversible controls 

6. induced velocities and flight velocity 

It is hoped that the present treatment, which is based upon 
linearized Euler' a equations, may suggest interesting alternatives 
to the more common hinge angle representation in other problems of 
rotor dynamics. 
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CONCLOBIONB 


1. For a helicopter "blade in which the restoring farces 
In vibration are due to Inertia effects rather than springs or 
structural stiffness, a simplified theory leads to conditions 
far the occurrence of flutter that can be represented, by a single 
generalized chart. For all examples In which the coning angle Is 
determined by a balance between lift forces and Inertia farces the 
critical condition corresponds to the same point on the chart. 

2 . The theory Indicates the dependence of flutter Instability 
upon the physical characteristics and the caning angle and pitch 
setting of the blade. 

3. The results show that the stability Is extremely sensitive 
to chordwlse position of center of mass with respect to aero- 
dynamic center. Forward movement of the center of mass Increases 
the stability. 

4. This paper 13 to be considered as only a preliminary study 
of blade flutter and as a starting point for further developments . 
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Figure 1. 
Inertia 


Principal Axes of 
of a Blade. 
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Figure 2 — Assumed Instantan- 
eous Force on Airfoil. 
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F/fure 3.- Stability Chart . 
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